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Abstract 

We report the Suzaku results of HESS J1614— 518, which is the brightest extended TeV gamma-ray 
O ■ source discovered in the Galactic plane survey conducted using the H.E.S.S. telescope. We discovered 

three X-ray objects in the field of view of the X-ray Imaging Spectrometer (XIS), which were designated 
as Suzaku J1614-5141 (src A), Suzaku J1614-5152 (src B), and Suzaku J1614-5148 (src C). Src A is an 
extended source located at the peak position of HESS J1614— 518, and therefore it is a plausible counterpart 
to HESS J1614— 518. The X-ray flux in the 2-10 keV band is 5xl0~ 13 erg s _1 cm -2 , which is an order of 
magnitude smaller than the TeV flux. The photon index is 1.7, which is smaller than the canonical value 
of synchrotron emissions from high-energy electrons found in some supernova remnants. These findings 
present a challenge to models in which the origin of the TeV emission is the inverse Compton scattering 
of the cosmic microwave background by accelerated electrons that emit X-rays via synchrotron emission. 
Src B is located at a relatively dim region in the TeV band image; however, its hydrogen column density 
is the same as that of src A. Therefore, src B may also be physically related to HESS J1614— 518. Src C is 
a foreground late-type B star. We also discovered a soft extended X-ray emission near HESS J1614— 518. 
Key words: acceleration of particles — X-rays: individual (HESS J1614— 518) — ISM: cosmic rays 
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1. Introduction At present, however, X-ray observations on these TeV 

sources are still limited. Matsumoto et al. (2007) ob- 
One of the striking results of the TeV 7-ray astron- served HESS J1616— 508 deeply with the Suzaku satel- 
omy with the H.E.S.S. Cerenkov Telescope is the discov- lite (Mitsuda et al. 2007), in addition to analyzing the 
ery of many very high-energy (VHE) objects along the archival data obtained from XMM-Newton. However, 
Galactic plane (Aharonian et al. 2005a; Aharonian et al. they found no X-ray counterpart with a severe upper 
2006). Most of these objects are spatially extended to limit of 3xl0 -13 erg s _1 cm -2 in the 2-10 keV band (see 
>0.°1, and could therefore be Galactic objects, such as su- also Landi et al. 2007). HESS J1303— 631 (Aharonian 
pernova remnants (SNR) or pulsar wind nebulae (PWN). et al. 2005b) is the same example. On the other hand, 
TeV 7-ray emissions are also detected from a young stel- HESS J1804— 216 has some X-ray counterparts (Landi 
lar cloud, Westerlund 2 (Aharonian et al. 2007). These et al. 2006; Cui & Konopelko 2006; Bamba et al. 2007; 
provide direct evidence that the inner Galactic region pro- Kargaltsev et al. 2007a); however, their physical connec- 
duces high-energy particles with energies above TeV level, tion to HESS J1804— 216 is still unclear. Recently, an ex- 
Observations of other wavelengths, particularly X-rays, tended X-ray emission from TeV J2032+4130 (Aharonian 
provide a better understanding of which particles are ac- et al. 2002) that has a peak position coincident with 
celerated and the relevant dynamics. High-energy elec- that of the 7-ray emission was discovered as a result 
trons predominantly emit synchrotron X-rays in a typical of a deep observation with XMM-Newton (Horns et al. 
interstellar magnetic field with strength of a few micro- 2007). However, Butt (2007) was skeptical about the de- 
Gauss. Therefore, the flux ratio of the X-rays to TeV tection of the diffuse emission, and argued that it can be 
7-rays provides valuable information that can distinguish explained by integrated emission of faint X-ray sources 
whether the accelerated particles are protons or electrons, detected with Chandra (Butt et al. 2006). Recently, an 
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extended X-ray emission from HESS J1834— 087 was dis- 
covered with XMM- Newton (Tian et al. 2007), which is 
a clear identification. However, there are also point-like 
sources in addition to the extended emission in the spec- 
tral region reported by Tian et al. (2007); therefore, the 
physical connection between the extended emission and 
HESS J1834— 087 remains unclear. Absolutely, it is im- 
perative to have an access to more X-ray samples with 
a physical connection to the VHE objects. Therefore, 
we observed HESS J1614-518 (hereafter HESS J1614) 
with the Suzaku satellite, which is the brightest among 
the newly discovered HESS objects (Aharonian et al. 
2005a; Aharonian et al. 2006). Two X-ray objects were 
found in the vicinity of HESS J1614 using the Swift /X- 
Ray Telescope (XRT) (Landi et al. 2006). However, nei- 
ther of them coincided with the 7-ray peak of HESS J1614; 
therefore their physical connection to HESS J1614 could 
not be clarified. Since HESS J 16 14 has a spatial exten- 
sion of ~ 10', high-sensitivity observations in a hard X- 
ray band with low and stable backgrounds are essential. 
The Suzaku X-ray Imaging Spectrometer (XIS; Koyama 
et al. 2007), combined with the X-Ray Telescope (XRT; 
Serlemitsos et al. 2007), satisfies these requirements and 
is regarded as one of the best instruments to search for 
diffuse X-ray emissions, particularly in the hard-energy 
band. Uncertainties quoted in this study are at the 90% 
confidence level, and errors on the data points in the X- 
ray spectra and radial profiles are at the la confidence 
level, unless otherwise stated. 

2. Observations and Data Reduction 

The brightest part of HESS J1614 was observed on 2006 
September 16. The on-source observation was followed by 
a background observation at an offset position (figure 1), 
where no known bright X-ray or 7-ray source existed in the 
field of view. The latitude of the offset position is selected 
to be almost the same as the position of HESS J1614, so 
that the Galactic ridge emission at the source position can 
be reliably subtracted (e.g., Worrall et al. 1982; Warwick 
et al. 1985; Koyama et al. 1986; Yamauchi, Koyama 1993; 
Kaneda et al. 1997; Ebisawa et al.. 2001; Ebisawa et al. 
2005; Sugizaki et al. 2001; Tanaka 2002; Revnivtsev et al. 
2006). The observations are summarized in table 1. 

The observations were performed with the four CCD 
cameras (XIS) located at the focal planes of four XRTs. 
One of the XIS sensors (XIS1) has a back- illuminated (BI) 
CCD, while the other three sensors (XIS0, 2, and 3) utilize 
front-illuminated (FI) CCDs. The XIS was operated in 
the normal clocking mode with no charge injection tech- 
nique. Although the non-imaging Hard X-ray Detector 
(HXD; Kokubun et al. 2007; Takahashi et al. 2007) was 
also available, in this study, we focus on the XIS data 
analysis because our main objective is searching for X-ray 
counterparts of HESS J1614 using the imaging capability. 

Data were reduced and analyzed using the processed 
data of version 1.2.2.3 1 , the HEADAS software version 

1 See http:/ /www. astro. isas.jaxa.jp/suzaku/process/history/vl223. 



6.2, and a calibration database (CALDB) released on 2007 
April 9. All data affected by the South Atlantic Anomaly 
and/or telemetry saturation were excluded. We also ex- 
cluded the data obtained with low elevation angles from 
the Earth rim that were below < 5° 2 . Finally, we removed 
the hot and flickering pixels. After these data screenings, 
the effective exposures were 44.5 ks and 53.1 ks on the 
HESS J1614 and the offset background, respectively. 

We checked the energy scale of each XIS sensor using the 
calibration source ( 55 Fe) and confirmed that the observed 
center energies of Ka and K/3 lines were in good agreement 
with the expected values (5.895 and and 6.490 keV) within 
the nominal calibration uncertainty of 0.2% (Koyama et 
al. 2007). 

3. Analysis and Results 

3.1. XIS image 

We extracted XIS images from each sensor using the 
screened data for the soft- and hard-energy bands. For 
the FI sensors, the soft- and hard-bands are defined as 
0.4-3.0 keV and 3.0-10.0 keV, respectively, while those 
for the BI sensor are defined as 0.3-3.0 and 3.0-7.0 keV, 
respectively. We excluded the corners of the CCD chips 
illuminated by the 5 Fe calibration sources. Images of 
the non- X-ray background (NXB) were obtained from the 
night Earth data provided by the XIS team (Tawa et al. 
2007) and subtracted from the HESS J1614 images. Then, 
the soft and hard images were divided by flat sky images 
simulated at 1.49 and 6.0 keV using the XRT+XIS simula- 
tor xissim (Ishisaki et al. 2007) for vignetting corrections. 
The images from the three FI sensors were summed and 
re-binned by a factor of 8. 

The XIS FI images of the HESS J1614 region shown 
in figure 2 were smoothed using a Gaussian function 
with a sigma of Of 28. The BI images were essentially 
the same except for the poorer statistics. In the hard- 
band image, an extended object with a peak position of 
(l,b) — (331.°64, — 0.°51) 3 was found, and was designated 
as Suzaku J1614— 5141 (src A). The position uncertainty 
was Of 8. Src A was not conspicuous in the soft-band 
image. Another bright X-ray object, found at (l,b) — 
(331P45, — 0.°60) 4 , was also bright in the soft X-ray band. 
The position uncertainty was Of 8. We designated this ob- 
ject as Suzaku J1614— 5152 (src B). In the soft-band im- 
age, an X-ray source was found at (/,&) = (331.°58, — 0.°62) 

2 We did not screen the data with the elevation angles from 
the bright Earth rim in order to maximize the statistics, al- 
though the standard data processing done at ISAS/JAXA or 
NASA/GSFC usually use this screening. The bright Earth may 
contaminate the soft-band spectra. In our case, however, we 
confirmed that the standard screening reduced the effective ex- 
posures but did not affect the XIS spectra. For example, exclud- 
ing the data with low elevation angles from the bright Earth rim 
of < 20° changed the counting rate of XIS1 in the 0.2—2 keV band 
from 0.1722 cts s _1 to 0.1726 cts s _1 , while the exposure time 
decreased from 44.5ks to 40.3ks. 

3 (a,<5)j2000.0 = (16 h 14 m 34 s , -51°41'00") 

4 (a,5)j2000.0 = (16 h 14 m 06 s ,-51°52'38") 
html 
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Table 1. Log of Suzaku observations. 



Name 


OBSID 


Pointing 
I 


direction 
b 


Observation start (UT) 


Effective exposure (ks) 


HESS J1614 

offset 


501042010 
501043010 


331.°5717 
330.°3961 


-0.°5274 
-0.°3760 


2006/09/15 16:00 
2006/09/16 11:02 


44.5 
53.1 




Fig. 1. Suzaku fields of view (thick boxes) overlaid on the H.E.S.S. smoothed excess map (same as figure 9 in Aharonian et al. 
2006). The scale bar below the figure represents the excess. The coordinates on the exterior frame are Galactic, while the grid shows 
the equatorial coordinates (J2000.0). The contour lines denote the radio emission (843 MHz) adopted from the Sydney University 
Molonglo Sky Survey (SUMSS; Bock et al. 1999). The contours are drawn from 0.1 Jy beam -1 on a logarithmic scale with a factor 
of 0.158, where the FWHM beam size is ~ 0.'7 x Of 9. Regions outlined by thin solid lines are the ones used in our spectral analyses 
(see section 3). 



5 and was designated as Suzaku J1614— 5148 (src C). The 
position uncertainty was 0'6. 

We obtained the radial profiles and compared with a 
point-spread function (PSF). As for the PSF, we ob- 
tained the radial profiles using the SS Cyg data observed 
on 2005 November 2 (OBSID=400006010), which is the 
verification phase data for the imaging capability of the 
XRT (Serlemitsos et al. 2007). Since the energy de- 
pendence of the PSF is almost negligible (Serlemitsos et 
al. 2007), it was extracted from the 0.4-10 keV band. 
Figure 3 shows the radial profile of src A in the hard- 
energy band. The profile cannot be fitted with the PSF 
plus a constant component model (x 2 /d.o.f. = 45.7/13), 
and therefore src A should be an extended source. Since 
src B is located near the edge of the XIS field of view, and 
therefore the PSF is largely distorted, it is impossible to 
study the radial profile of src B in detail. However, the 
hard-band image suggested that src B is extended. The 
radial profile of src C is consistent with the PSF plus a 

5 (a,5)j20oo.o = (16 h 14 m 47 s ,-51°48 / 36 // ) 



constant component (x 2 /d.o.f. = 37.1/26), supporting the 
assumption of a point source. 

3.2. XIS spectrum 

3.2.1. Background spectra: blank and offset positions 

Since the source position is on the Galactic ridge, the 
local Galactic emission (the ridge emission) , as well as the 
NXB and cosmic X-ray background, cannot be ignored. 
To account for the background, we compared the spec- 
trum from the blank region near the sources (see the solid 
line in figure 2) with that from the offset position. The 
spectra of the offset position correspond to a circular re- 
gion of 6' radius at the center of the field of view. 

For the most accurate NXB estimation, we sorted the 
night Earth data so that the cut-off rigidity distribution 
was the same as that of the HESS J1614 and offset ob- 
servations. Since the NXB depends on a location on the 
detector (Yamaguchi et al. 2006), we extracted the NXB 
spectra from the same regions used to determine the blank 
and offset regions in detector coordinates (DETX/Y). The 
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Fig. 2. Suzaku XIS FI (XISO+2+3) images of the HESS J1614 region in the Galactic coordinates: (a) 0.4-3 keV and (b) 3-10 keV 
band. The images were smoothed using a Gaussian function with a sigma of 0.'28. Vignetting correction was applied after subtracting 
NXB as described in the text. Regions represented with solid lines are used in our spectral analyses (see section 3.2). Plus marks 
represent the peak position of each source. 
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Fig. 3. Radial profile of src A extracted from the hard-band 
image of the XIS FI sensor (XISO+2+3). The solid line rep- 
resents the XIS PSF profile with a constant component. The 
downward arrows represent the positions of src A2 and A3 
defined in section 3.2.3 (see figure 7). 

NXB spectra thus obtained were then subtracted, and the 
resulting spectra of the blank and offset regions are shown 
in figure 4. The offset spectra were normalized to the same 
area as that of the blank region. We can clearly observe 
that the spectra of the blank region show an enhancement 
less than 3.0 keV. 

The excess spectra of the blank region after subtracting 
the data of the offset position are shown in figure 5. We 
fitted these spectra with an absorbed power-law model, 
and the cross sections of photoelectric absorption were 
obtained from Morrison & McCammon (1983). We ob- 



tained detector responses (RMF) and telescope responses 
(ARF) using both the xisrmf gen and xissimarf gen soft- 
ware (Ishisaki et al. 2007). In this hydrogen- 
equivalent column density (A%), a photon index (T), and 
normalization are free parameters. The best-fit parame- 
ters are summarized in table 2. Although the data were 
fitted acceptably with the power-law model, we can see 
some indications of emission lines in the residual, par- 
ticularly at ~1.3 keV. When we added a Gaussian line 
model with a — (i.e., the delta function) to the power- 
law model, the peak energy became 1.331q'o4 keV and the 
equivalent width was 24.0 (<53.1) eV, although the F-test 
does not require the use of a Gaussian model. These re- 
sults suggest the emission line of MgXI and a thermal 
spectrum. Therefore we tried to fit the spectra with a 
thin-thermal plasma model (the APEC model; Smith et 
al. 2001), where the solar abundance of elements is ob- 
tained from Anders, Grevesse (1989). The free parameters 
are plasma temperature (kT), metal abundance, A%, and 
normalization. The spectra were fitted acceptably with 
the plasma model, as summarized in table 2. However, 
the metal abundance is almost zero. 

Since the blank region is in a dim area of both the soft- 
and hard-energy bands (figure 2), the soft excess must 
also exist in the regions of src A, B, and C. As for the 
source spectra, we concern the excess flux over this soft 
excess emission, and therefore used the blank region as 
background in the following analyses of the source spectra. 
3.2.2. Extraction of the source spectra 

We represented the source regions of src A, B, and C 
as solid lines in figure 2. We extracted light curves corre- 
sponding to the sources in the soft- and hard-energy bands 
from these regions and observed that none of them exhib- 
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Fig. 4. XIS spectra of the blank and offset regions: (a) FI spectra (XISO+2+3) and (b) BI (XIS1) spectra. The offset spectra were 
normalized to the same area as that of the blank region. The NXBs were subtracted as described in the text. 
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Fig. 5. Spectra of the blank region and the best-fit power-law 
model. The spectra of the offset region are subtracted as back- 
ground. Black and gray lines represent the data and model 
for the XIS FI (XISO+2+3) and BI (XIS1), respectively. 

ited clear time variability. We extracted the XIS spectra 
of src A, B, and C by subtracting the NXB spectra, in the 
same manner as in the previous section. We then com- 
bined the NXB subtracted spectra obtained from three 
FI sensors. 
3.2.3. Src A 

The spectra of src A are shown in figure 6. Although 
the statistics are limited, no clear emission line was ob- 
served. We therefore fitted the spectra with the absorbed 
power-law model and obtained acceptable results. The 
best- fit parameters are summarized in table 2. In ad- 
dition, we also tried to fit the spectra with the ther- 
mal plasma model. The best-fit parameters were ATh = 
1.04±g;|Jxl0 22 cm" 2 , kT = 10.0l 26 2 ° keV, and an abun- 
dance of 0.18(< 0.93) solar with x 2 /d.o.f. = 80.01/67. 
From a statistical viewpoint, neither of the models is re- 
jected at a confidence level of 90 %. However, the thermal 
plasma model yields an uncomfortably high temperature 



Fig. 6. XIS spectra of src A and the best- fit power-law 
model. Black and gray lines represent the data and the model 
for the XIS FI (XISO+2+3) and BI (XIS1), respectively. 

and an extremely low abundance, and therefore is not 
practicable. We therefore adopted the power-law model 
in all subsequent analyses and discussions. 

Figure 7(a) is a close-up view of src A in the XIS 
FI hard-band image (figure 2). One could argue that 
src A is not a single extended source, but instead con- 
sists of multiple objects. To verify this possibility, we 
defined three circular regions designated as Al, A2, and 
A3 in the figure, and extracted the corresponding spec- 
tra. The radii of the source regions are 2', 1', and V 
for Al, A2, and A3, respectively. However, the statis- 
tics are poor and with the exception of the corresponding 
normalizations, no clear difference between these sources 
was observed (figure 7(b)). The spectra are described by 
a power- law model with a best-fit x 2 /d.o.f. of 25.72/31, 
13.35/23, and 23.03/23, for Al, A2, and A3, respec- 
tively. The observed fluxes in the 0.5-10.0 keV band are 
2.9xl0~ 13 erg s -1 cm -2 , l.lxlO -13 erg s -1 cm -2 , and 
0.95xl0~ 13 ergs -1 cm -2 , for Al, A2, and A3, respec- 
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Fig. 7. (a) Close-up view of src A in the XIS FI hard-band 
image (figure 2). A large ellipse shows the spectral region for 
src A. Dotted circles designated as Al, A2, and A3 are the 
regions we studied to check for spatial variation of the spectra 
in src A. (b) Spectra extracted from regions Al, A2, and A3, 
as shown in figure 7(a). For clarity, only the spectra obtained 
by the FI sensors (XISO+2+3) are shown. 



tively. The best-fit parameters are summarized as fol- 
lows: N u = 1.40tg;!?x 10 22 cm" 2 and T = 1.41±g g for the 
Al region, N u = 0.66t°;;gx 10 22 cm" 2 and V = 1.39±g g 
for the A2 region, and A^ H = 1.47±o.7o><10 22 cm" 2 and 
r = 2.26^o - 58 f° r the A3 region. All these parameters are 
consistent with the best-fit values of src A presented in 
table 2, which supports the hypothesis that src Al, A2, 
and A3 are substructures of src A. 
3.2.4. Src B 

The spectra of src B are shown in figure 8. We fit- 
ted the power-law model with the best-fit parameters as 
listed in table 2. Although no emission line was ob- 
served in the spectra, we also tried fitting a thermal 
plasma model. The best-fit x 2 value of 196.13 for 159 de- 
grees of freedom rejected the validity of a thermal plasma 




Fig. 8. XIS spectra of src B, shown with the best-fit pow- 
er-law model. Black and gray lines represent the data and the 
model for the XIS FI (XISO+2+3) and BI (XIS1), respectively 



model at a confidence level of 97%. Nevertheless, for the 
sake of comparison, following are the best-fit parameters: 

N n = 0.84i2;5i xl ° 22 cm " 2 > kT = L28 -o:i3 keV > and a 
metal abundance of zero (< 0.019) solar. 

Src B is at the edge of the field of view. According 
to a simulation performed using xissim, the probability 
of detecting incident photons reduces by a factor of ~0.7 
because some photons fall outside of the field of view. 
Furthermore, considering the fluctuation of the pointing 
direction due to the variation of the relative alignment 
between the XRT and the Attitude and Orbit Controlling 
System (AOCS; Serlemitsos et al. 2007), even more pho- 
tons can escape from the field of view. Consequently, 
all the obtained parameters, particularly flux, could have 
large systematic errors. 
3.2.5. Src C 

We found that the spectra of src C can be fitted with 
the thermal plasma model, as shown in figure 9. The best- 
fit parameters are listed in table 2. Although the spectra 
can also be fitted with a power-law model having best-fit 



1.01 



+0.53. 
-0.36 



xlO 22 cm 2 and V 



8.635;1 



values as N-r 

with x 2 /d.o.f. = 63.76/52, the obtained photon index is 
unreasonably large. Therefore, the thermal plasma model 
seems more suitable for practical use. 

4. Discussion 

4.1. Src A 

Src A is the closest source to the peak position 
of HESS J1614, and both the X-ray and TeV 7-ray 
emissions are extended. Therefore, src A is a plau- 
sible X-ray counterpart to HESS J1614. The X-ray 
spectrum is described using the power-law model with 
r = 1.73, and the absorption-corrected flux in the 2- 
10 keV band is 5.3x 10 _13 erg s _1 cm -2 , while the 7- 
ray spectrum yields T=2.46 and the flux in the 1— 
10 TeV band from a circular region of 0.°4 is estimated 
to be 1.8xl0 -11 ergs -1 cm -2 (Aharonian et al. 2006). 
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Table 2. Best-fit results of the XIS spectra. 
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Model* 


PL 


APEC 
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PL 


"DT 

PL 


APEC 


at /1 n22 — 2\ 

iVH (10 cm ) 


24+ - 08 


10+ 005 

u - lu -0.04 


i r.1+0.50 

1.21I .4i 


1 /i +0.14 


n OO+0.19 

0.281 0>ig 


T/kT (keV) 


Q OO+0.19 

z -°°-0.17 


9 79+O.5O 
z - ' z -0.43 


-, 70+0.33 

'°-0.30 


6 n+0.24 
O.DU_ 22 


0.63^ 


Abundance (solar) 




0.00(<0.029) 






O.i8± ;lg 


pobs If 
^O.S-lOkeV 


0.10 


0.096 


5.4 


5.1 


0.86 


pabscor § 
^O.S-lOkeV 


0.15 


0.11 


8.4 


35.8 


2.0 


pobs f 
- r 2-10keV 


0.058 


0.050 


4.8 


3.0 


0.071 


pabscor § 
" r 2-10keV 

X 2 /d.o.f. 


0.059 


0.051 


5.3 


3.6 


0.076 


121.45/103 


111.37/102 


78.38/68 


172.65/160 


42.45/51 



f Background spectrum used for the analysis. 

\ Model used for the spectrum fitting: "PL" is a power-law model, and "APEC" is a thin thermal plasma model. 
% Observed flux in the 0.5-10 keV or 2-10keV band in the unit of 10 -13 erg s _1 cm -2 . For the blank region, values 
are normalized to a region of 1 arcmin 2 (i.e., erg s _1 cm -2 arcmin -2 ). 

§ Absorption corrected flux in the 0.5-10 keV or 2-10keV band in the unit of 10 -13 erg s _1 cm -2 . For the blank 
region, values are normalized to a region of 1 arcmin 2 (i.e., erg s _1 cm -2 arcmin -2 ). 
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Fig. 9. XIS spectra of src C, shown with the best-fit thermal 
plasma model. Black and gray lines represent the data and the 
model for the XIS FI (XISO+2+3) and BI (XIS1), respectively 

Therefore, the flux ratio F(l - 10 TeV)/F(2 - 10 keV) is 
~34, which is one of the largest values observed among 
extended VHE objects (see Matsumoto et al. 2007 and 
references therein). 

We first examined the difference between the peak posi- 
tions of src A and HESS J1614. Since the peak position of 
TeV 7-ray is not presented in the study by Aharonian et 
al. (2006), we only worked in the bright region of figure 1 
and fitted it with a Gaussian function. The peak position 
thus estimated is (l,b) = (331.°65, -0.°52) 6 , which is ~0.'8 
away from the X-ray peak of src A. Considering that the 
best H.E.S.S. angular resolution is ~5' (Aharonian et al. 
2006), and from the position uncertainty of src A (~0.'8), 
we can conclude that src A spatially coincides with the 

6 (a,5)j 20 oo.o = (16 h 14 m 41 s , -51°41'03") 



peak of HESS J1614. 

If src A is physically associated with HESS J1614, the 
best-fit column density (1.2xl0 22 cm -2 ) suggests that the 
distance to HESS J1614 is roughly ~10 kpc, since the 
best-fit value is approximately equal to half of the to- 
tal Galactic HI column density towards the HESS J1614 
region (~ 2.2 xlO 22 cm" 2 ; Dickey & Lockman 1990). 
Assuming a distance of 10 kpc, the luminosity of src A 
in the 2-10 keV band is 6xl0 33 erg s -1 . 

This region was also observed with the Swift 
XRT (Landi et al. 2006); however, no source was found 
at the position of src A. This is probably due to the lim- 
ited exposure time (~ 1700 s), and/or the small effective 
area of the Swift XRT. Furthermore, src A is near the 
edge of the field of view of the Swift XRT, which makes 
the detection more difficult. No counterparts in the other 
wavelengths are found in literatures. 

If the origin of the TeV emission is the inverse Compton 
scattering of the cosmic microwave background by high- 
energy electrons of a single population, the large flux ratio 
F(l - 10 TeV) /F (2 -10 keV) ~34 requires a low magnetic 
field. The spectral energy distribution (SED) of src A and 
HESS J1614 are plotted in figure 10. In the same fig- 
ure, we also plotted the estimated synchrotron flux of the 
high-energy electrons in the magnetic fields with intensi- 
ties B =0.1, 1, and 10 //Gauss. The low X-ray flux of src A 
requires a magnetic field of B ^ 1//G, which is smaller 
than the typical interstellar magnetic field intensity of a 
few micro- Gauss. This is similar to what is observed in 
HESS J1616-508 (Matsumoto et al. 2007); however, but 
in that possibility of a magnetic field of a few 

micro- Gauss and a strong cut-off in the electron energy 
distribution greater than ~ 10 14 eV cannot be excluded. 
However, Src A has a small spectral index, which does not 
favor this possibility. 

On the other hand, interstellar radiations whose wave- 
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lengths are shorter than those of the CMB (e.g., far- 
infrared) also have large energy density (Moskalenko et 
al. 2006). If the seeds of the inverse Compton are such 
photons, the electrons produce synchrotron emission at 
wavelengths longer than X-rays in regions with a mag- 
netic field of a few micro-Gauss, and therefore our Suzaku 
results would have no constraints on the magnetic field 
intensity. 

Since src A is extended, it may be an SNR. However, 
the observed X-rays are probably not due to ordinary 
synchrotron emission, because the photon index (~1.7) is 
smaller than the typical values of non-thermal SNRs such 
as SN1006 (~2.5-3.0; Bamba et al. 2005). The X-ray spec- 
trum of the SNR 7 Cygni also has a small photon index of 
0.8-1.5 (Uchiyama et al. 2002), which can be considered 
as a non-thermal bremsstrahlung from accelerated elec- 
trons; the energy distribution of the accelerated electrons 
below the Coulomb break is flattened due to the Coulomb 
interactions with ambient molecular clouds (Uchiyama et 
al. 2002). In that scenario, a photon index less than 1.5 
is expected for viable acceleration mechanisms. Since the 
photon index of src A is marginally consistent with this 
criterion, the observations of src A could be explained with 
this scenario. In addition, undetected molecular clouds re- 
sponsible for the bremsstrahlung may exist. Yamazaki et 
al. (2006) suggested that the TeV-to- X-ray flux ratio of 
old SNRs with an age of ~ 10 5 years can be very large, 
in some cases more than 100. The small photon index of 
src A may support this scenario, where X-ray emission is 
dominated by synchrotron emission from secondary elec- 
trons that may have a smaller index than that of the pri- 
mary electrons. However, no radio counterpart has been 
observed in the HESS J1614 region, which could pose a 
severe problem for the SNR scenarios. For example, there 
is no enhancement in the 843 MHz band (figure 1), where 
the rms noise level is ~2 mJy arcmin -2 (Bock et al. 1999). 

Massive stars can be the origin of the extended TeV 7- 
ray emission; Torres & Domingo- Santamaria (2007) pro- 
posed that stellar winds can produce hadronic 7-rays by 
interacting with cosmic rays. Anchordoqui et al. (2007) 
argue that TeV 7-rays are the result of Lorentz-boosted 
MeV 7-rays, which are emitted due to the de-excitation 
of daughter nuclei resulting from collisions between high- 
energy nuclei (the cosmic rays) and UV photons from the 
massive stars. In either case, X-rays may originate from 
stellar winds from the OB stars. These models have been 
proposed in order to explain the origin of the unidenti- 
fied TeV source, TeV J2032+4130. Although the X-ray 
photon index of TeV J2032+4130 is similar to that of 
src A, the flux ratio F(l - 10 TeV)/F(2 - 10 keV) (~ 
2-3) is much smaller than that of src A (Horns et al. 
2007). Furthermore, no OB stars have been observed 
around HESS1614. 

Another possibility is a PWN; some PWNe have a small 
photon index (Fleishman & Bietenholz 2007 and refer- 
ences therein). The luminosity, assuming the distance 
of 10 kpc, is not unusual for X-ray PWNe (Cheng et 
al. 2004). A pulsar driving the nebula may be hidden 
in the middle of src A, or src B itself could be the pul- 
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Fig. 10. Spectral energy distribution of src A and 
HESS J1614 from the X-ray to TeV 7-ray bands. The syn- 
chrotron radiation from accelerated electrons, which kick the 
3K background up to the TeV energy range, is plotted toward 
the left for three different values of the magnetic field inten- 
sity. The regions used to estimate the SED in the TeV 7-ray 
and X-ray bands are different. 

sar. However, the flux ratio F(l - 10 TeV)/F(2 - 10 keV) 
is unusually large compared with the typical PWNe; 
the ratios of Crab and RCW89 are 2.6 xlO -3 and 0.41, 
respectively (Willingale et al. 2001; Aharonian et al. 
2004; Aharonian et al. 2005c; DeLaney et al. 2006). In 
addition, the lack of radio emissions may be difficult to 
explain by assuming a PWN origin. 

4.2. Src B 

The best-fit value of the column density suggests that 
the distance to src B from us is similar to that to src A. 
There is an optical (Monet et al. 1999) and an infrared 
counterpart (Skrutskie et al. 2006), 2MASS J16140610- 
5152264, at (l,b) = (331.°4579,-00.°5973), whose color also 
suggests extinction (Landi et al. 2006). Assuming a dis- 
tance of 10 kpc, the luminosity of src B in the 2-10 keV 
band is 4xl0 33 erg s _1 . 

Src B was also detected with the Swift XRT (Landi 
et al. 2006). Their observed count rate was (6.61 ± 
2.11) xlO -3 c s" 1 in the 0.2-10 keV band. Landi et 
al. (2006) converted the count rate to the absorption- 
corrected flux of 5.4x 10 _13 erg s _1 cm -2 in the 2- 
10 keV band by assuming a Galactic column density of 
2xl0 22 cm -2 (Dickey & Lockman 1990) and a Crab- like 
spectrum. However, these parameters are not supported 
by our Suzaku results. We converted the count rate to the 
flux in the 2-10 keV band by assuming the best-fit param- 
eters listed in table 2 with "WebPIMMS" 7 . The resulting 
observed flux was (1.43±0. 46) x 10 -13 erg s _1 cm -2 . This 
suggests that time variability exists in src B, although the 
Suzaku flux of src B may have uncertainty, as described 
in section 3.2. 

Since src B is in a relatively dim region in the TeV 7- 

7 http: / /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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ray band, this object seems hardly related to HESS J1614. 
However, if src A is the X-ray counterpart of HESS J1614, 
there still may be a remote possibility that src B may 
also be physically connected to HESS J1614, as suggested 
by the best-fit column density. One may speculate that 
HESS J1614 is an asymmetric PWN, such as MSH 15- 
52 (Aharonian et al. 2005c), and src B is a pulsar driv- 
ing the nebula. The TeV 7-ray morphology suggests that 
HESS J1614 may consist of double-sided structures, and 
src B seems to reside between them (figure 1). The im- 
age of the TeV 7-ray emission from MSH 15-52 is also 
elongated and there is an indication that the TeV emis- 
sion has peaks on both sides of pulsar PSR B 1509— 58. 
Pulsars within nebulae having a luminosity of 10 33 erg s _1 
are not unusual. The spin down luminosity should be 
~ 10 37 erg s" 1 (Cheng et al. 2004), then the TeV 7-ray 
emission of HESS J1614 can be driven by a few percent 
of the spin-down luminosity. The extremely large pho- 
ton index of 3.6 is enigmatic, but consistent with the 
X-ray photon index of PSR J1809-1917 (Kargaltsev & 
Pavlov 2007b), which could conceivably associated with 
HESS J 1809— 193. However, at present, we observe no 
pulsation or radio emission from src B, as shown in fig- 
ure 1. 

4.3. Src C 

The best-fit column density suggested that src C is 
a foreground source. From the SIMBAD Astronomical 
Database operated at CDS, Strasbourg, France 8 , we found 
a counterpart at (l,b) = (331.°5819, -00.°6251) 9 , which is 
a B9V star: HD 145703. The best-fit temperature and 
abundance are consistent with the X-ray observations of 
other late- type B stars (Stelzer et al. 2003). 

4-4- Soft emission in the blank region 

Within the bounds of our statistics, it is unclear 
whether the emission observed in the blank region is ther- 
mal or non-thermal. Given its soft and unabsorbed na- 
ture, this source may not be connected to HESS J1614 at 
all. Instead, we suspect that the emission is due to spatial 
fluctuation of the Galactic ridge emission, which is com- 
prised of high (~ 7 keV) and low (~ 0.7 keV) temperature 
components (Kaneda et al. 1997). Since the best-fit col- 
umn density is small, the soft excess emission may include 
the low temperature component. The best-fit temperature 
is between the high and low temperature components, and 
therefore the soft excess may also include high tempera- 
ture component. We found no evidence suggesting the 
presence of soft emission in the hard band (0.5-2.0 keV) 
image of the ROSAT All-Sky X-ray Survey Broad Band 
data (RASS3bb, RASS-Cnt Broad). 

5. Conclusions 

The Suzaku results on the VHE 7-ray object 
HESS J1614— 518 are summarized below. 

8 http://simbad.u-strasbg.fr/simbad 

9 («,5)j2ooo.o = (16 h 14 m 48F5 1 -51°48 / 31."3) 



1. Three X-ray objects, Suzaku J1614— 5141 
(src A), Suzaku J1614-5152 (src B), and 
Suzaku J1614— 5148 (src C) are found in the 
Suzaku XIS field of view. 

2. A hard extended source, src A, is the best candidate 
for the X-ray counterpart of HESS J1614— 518. 

3. Src A has a large flux ratio of F(l - 10 TeV) /F(2 - 
10 keV) ~ 34 and a relatively small X-ray photon 
index of 1.7. It is difficult to explain these fea- 
tures by assuming that the TeV emission is due to 
the inverse-Compton scattering of the cosmic mi- 
crowave background and that the X-rays are due to 
synchrotron radiation emitted by high-energy elec- 
trons. 

4. The small photon index of src A can be explained in 
terms of non-thermal bremsstrahlung emission from 
the loss-flattened electron distribution. 

5. Src B may not be physically related to 
HESS J1614— 518. However, one remote possi- 
bility is a PWN scenario, in which src B is a pulsar 
driving the nebula. 

6. Src C is a foreground B9V star, HD 145703. 

7. We found a soft excess emission near 
HESS J1614— 518, which is characterized ei- 
ther by a power- law with r = 2.4 or a thermal 
plasma with kT — 2.7 keV. This is probably not 
directly related to HESS J1614-518, but would 
nonetheless form a part of the Galactic ridge 
emission. 
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